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SUMMARY 


Several  methods  of  obtaining  the  time  response  of  linear  systems  to 
either  a  unit  Impulse  or  an  arbitrary  input  from  frequency-response  data 
are  described  and  compared. 

Comparisons  indicate  that  all  the  methods  give  good  accuracy  when 
applied  to  a  second-order  system;  the  main  difference  is  the  required 
computing  time.  Several  of  the  methods  when  applied  to  higher  order 
systems  require  excessive  computing  time  in  order  to  obtain  the  same 
degree  of  accuracy.  The  methods  generally  classified  as  inverse  Laplace 
transform  methods  were  found  to  be  most  effective  in  determining  the 
response  to  a  unit  impulse  from  frequency-response  data  of  higher  order 
systems . 

Some  discussion  and  examples  are  given  of  the  use  of  the  methods 
as  flight-data-analysis  techniques  in  predicting  loads  and  motions  of  a 
flexible  aircraft  on  the  basis  of  simple  calculations  when  the  aircraft 
frequency  response  is  known. 


INTRODUCTION 


The  frequency-response  type  of  analysis  used  on  linear  systems  has 
found  extensive  application  in  the  field  of  aircraft  stability  and  in 
the  determination  of  overall  dynamic  characteristics  of  an  aircraft.  In 
fact  the  current  trend  to  perform  analysis  on  flight  test  data  in  frequency- 
response  form  appears  to  be  gaining  favor  especially  in  the  case  of  flex¬ 
ible  aircraft.  Considerable  emphasis  has,  therefore,  been  given  to  methods 
of  determining  the  frequency  response  of  a  system  from  transient  responses. 
A  concise  re'sume  and  comparison  of  methods  for  obtaining  the  frequency 
response  from  transient  responses  are  presented  in  reference  1. 

This  paper  briefly  considers  the  next  step  in  the  process  -  the 
converting  of  the  information  contained  in  a  frequency  response  to  the 
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time  plane  in  the  form  of  the  response  to  a  unit  impulse.  For  certain 
purposes  this  conversion  gives  data  in  more  useful  form.  The  response 
of  a  linear  system  to  a  unit  Impulse  may  he  used  in  conjunction  with  * 

Duhamel's  (superposition)  integral  to  determine  (l)  the  aircraft  transient 
response  to  any  type  of  input  or  (2)  the  input  required  to  cause  any 
required  aircraft  transient  response. 

It  appears  that  methods  of  converting  frequency-response  data  to 
transient  data  present  a  flight-data-analysis  technique  which  permits 
the  prediction  of  aircraft  motions  and  loads  for  a  flexible  aircraft 
without  knowledge  of  the  equations  of  motion  relating  the  input  and 
output.  These  methods  also  bypass  the  need  for  computing  transfer- 
function  coefficients  or  stability  derivatives  in  predicting  these  loads 
and  motions.  Such  predictions  are  important  in  anticipating  the  motions 
and  loads  for  more  hazardous  aircraft  maneuvers.  These  methods  are  also 
useful  in  predicting  time  responses  of  complicated  linear  systems  whose 
frequency  response  is  known. 

The  purpose  of  this  report  is  to  collect  and  briefly  compare  a  few 
of  the  methods  now  available  for  performing  this  operation.  The  methods 
are  compared  on  the  basis  of  accuracy,  computing  time  required,  and 
applicability  of  the  method  to  higher  order  systems.  Some  extensions  of  4 

these  methods  are  also  given.  With  such  information  available,  engineers 
may 'then  select  the  method  which  best  fits  their  needs. 

V 


SYMBOLS 


F(t) 

H(ico) 

H(s) 

h(t) 


K,L 

M 


function  of  time 

frequency  response  of  a  system 

transfer  function  (in  terms  of  Laplace  variable  s) 
relating  input  and  output 

time  response  of  a  system  to  a  unit  impulse 

index  of  summation 

transfer  coefficients  defined  by  equation  (8) 

limits  on  summation  of  P-transform 
Mach  number 
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m 

N 


n 

P 

Re  [H(ico)] 
r 

s 

t 

X 

x(s) 

x(t) 

X 

5 

A 

5 

6(t) 

3 


T 


0 


03 


index  of  summation  of  P-transform 

total  number  of  data  points  used  in  solution  of 
equation  (9) 

limit  of  summation 

polynominal  transform  operator 

real  part  of  frequency  response  defined  by  equation  (21) 
ordinate  of  pulse  used  to  fit  Re|j3(ia))] 

Laplace  variable,  c  +  ico 

time,  sec 

output 

Laplace  transform  of  x(t) 

response  of  a  system  to  an  arbitrary  input 

amplitude  ratio  of  frequency  response 

indicates  increment 
input 

time  history  of  input  to  a  linear  system 
pitching  velocity,  radians/sec 
value  of  time,  sec 
phase  angle,  deg 

phase  angle  between  output  x  and  input  5  of  frequency 
response;  negative  phase  angles  indicate  lag 

circular  frequency,  radians/sec 

damped  natural  frequency,  radians/sec 

fundamental  frequency,  radians/sec 
undamped  natural  frequency,  radians/sec 
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Subscripts: 
e  elevator 

i  index  of  summation  indicating  row 

j  index  of  summation  Indicating  column 

^  A  tilde  indicates  polynominal  transform  of  function;  for  example, 
F(x)  denotes  polynominal  transform  of  F(t). 

DESCRIPTION  AND  DISCUSSION  OF  METHODS 


In  this  section  the  methods  for  determining  the  response  to  a  unit 
impulse  and  the  methods  for  obtaining  the  response  to  arbitrary  inputs 
are  discussed* 

The  methods  for  obtaining  the  time  response  to  a  unit  impulse  from 
frequency-response  data  fall  into  several  basic  categories  which  may  be • 
generally  classified  as  follows:  „ 

(1)  Inverse  Laplace  transform  methods 

(2)  Fourier  method  V 

(3)  Other  methods 

These  methods  are  briefly  outlined  to  indicate  the  techniques  involved. 

In  order  to  describe  the  computations  required  and  the  accuracies 
obtained,  each  method  is  applied  to  the  frequency-response  data  of  fig¬ 
ure  1  which  defines  a  simple  second-order  system  described  by  the 
transfer -function 


1 

s2  +  6s  +  10 


(1) 


It  is  usually  intended  that  these  methods  be  applied  to  higher  order 
systems  as  will  be  shown  later  in  the  paper.  No  attempt  is  made  to 
repeat  the  development  of  the  methods  since  this  information  may  be 
obtained  from  the  references. 


Inverse  Laplace  Transform  Methods 

Floyd  *  s  method . -  The  method  developed  by  George  F.  Floyd  and  described 
in  detail  in  reference  2  is  referred  to  as  Floyd’s  method.  Floyd  shows  * 
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that  the  inverse  Laplace  transform  h(t)  of  H(s)  given  by  the 
integral 


PC+ioo 

h(t)  =  ^—7  /  H(s)e  S  ds 


2rti 


(2) 


C-ioo 


is  for  all  positive  values  of  time  equivalent  to 


Re  [k(ia))]  cos  ten  du> 

0 


(3) 


This  operation  is  based  on  the  assumptions  that  H(s)  may  be  written 
as  the  ratio  of  two  rational  polynomials  in  s  with  real  and  constant 
coefficients,  that  lim  H(s)  =  0,  and  that  H(s)  has  no  poles  in  the 

s— y«> 

right  half  of  the  s -plane  or  on  the  imaginary  axis.  The  procedure 
for  performing  the  integration  required  by  equation  (3)  is  to  plot 
Re[H(ia))]  against  co  and  then  to  approximate  the  exact  shape  with  a 
series  of  straight-line  segments.  The  straight-line  approximation  is 
written  as  a  sum  of  trapezoidal  functions  and  equation  (3)  is  applied 
to  each  of  the  trapezoids;  the  resulting  time  functions  due  to  each 
trapezoid  are  then  added  to  obtain  h(t). 

A  simple  illustration  is  shown  in  sketch  A: 


Re 


The  time  function  associated  with  a  typical  trapezoid  (as  shown  by  the 
cross-hatched  section)  may  be  expressed  as 
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* 

where 


0*^ 

2  2 

%  -  % 

a2 - — 

A2  =  r2ad2 


In  general  for  n  trapezoids  the  time  response  is  given  by 


n 


h(t)  =  A. 


bin  oo^t' 

_ _ 

oo. t 
J 


sin  A-jt 

V 


(5) 


In  figure  2  the  real  part  Re[H(ia))]  of  the  simple  second-order 
system  computed  from  figure  1  is  shown  plotted  against  co.  The  func¬ 
tion  h(t)  is  then  easily  evaluated  since  tables  of  s:^  -  are  given 
in  reference  2.  In  the  upper  part  of  figure  2,  two  straight-line  fits  y 

to  the  Re  [H ( ico)]  are  shown,  one  a  five-line  fit  and  the  other  an  eleven¬ 
line  fit.  For  clarity,  however,  only  the  points  and  not  the  connecting 
lines  of  the  fit  are  shown.  In  the  lower  part  of  figure  2  the  resulting 
responses  to  a  unit  impulse  computed  by  Floyd 1 s  method  are  compared  with 
the  exact  response  to  a  unit  impulse.  The  accuracy  of  Floyd's  method 
depends  on  the  number  of  lines  used  to  fit  the  Re(j5(icD)3  and  the  loca¬ 
tion  of  the  cutoff  frequency.  The  cutoff  frequency  is  defined  as  the 
maximum  frequency  at  which  the  Re  [H(ioj)J  was  fitted.  It  should  be 
noted  that  the  limits  of  the  definite  integral  of  equation  (3)  are  0 
to  co;  however,  in  the  practical  case  the  Re[H(iaoJ|  is  cut  off  at 
some  finite  frequency.  This  error  is  reflected  at  the  low  values  of 
time,  especially  at  t  =  0.  For  a  given  number  of  lines,  the  accuracy 
also  depends  on  the  judicious  fit  of  the  lines. 

Numerical-integration  method.-  In  order  to  use  automatic  computing 
machines  to  perform  the  inverse  Laplace  transform  method,  the  necessary 
operations  indicated  by  equation  (3)  are  performed  in  the  following 
manner  by  using  a  numerical-integration  method.  For  a  given  value  of 
time  t,  the  product  curve  Re  [H(io>)]  cos  tea  is  evaluated  over  the  o> 
range  and  integrated  by  numerical-integration  techniques  which  give  one 
point  on  the  time  history  of  the  response  to  a  unit  impulse.  ISy  repeating 
the  above  computation  for  all  the  desired  values  of  time,  the  time  response  w 

can  be  obtained. 
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This  method  was  applied  to  the  frequency-response  data  of  figure  1. 
The  accuracy  of  the  method  is  shown  in  figure  3.  In  the  upper  part  of 
figure  3  the  error  due  to  the  interval  Ad  chosen  for  the  computation 
is  shown  by  the  comparison  between  the  circle  and  square  symbols.  The 
accuracy  of  the  computation  for  the  case  of  Ad  =1.0,  as  indicated  by 
the  circle  symbols,  was  not  satisfactory  beyond  t  =  2.0  seconds.  This 
result  is  to  be  expected  since  the  interval  Ad  =  1.0  was  too  large  to 
permit  numerical- integration  methods  to  perform  adequately  the  integra¬ 
tion  required  by  equation  (3). 

Rectangular -pulse  method .  -  A  method  for  determining  the  time  response 
to  a  unit  impulse  from  frequency -response  data  has  recently  been  given 
in  reference  3.  It  is  referred  to  herein  as  the  rectangular -pulse  method. 
The  method  involves  the  use  of  tables  for  time-plane  values  equivalent 
to  unit  rectangular  pulses  of  the  Re[k(ia)j]  .  The  method  requires  fitting 
the  Re[H(ia))]  with  a  series  of  rectangular  pulses  or  a  staircase  func¬ 
tion  so  that  the  area  under  the  curve  is  equal  to  the  area  of  the  pulse 
in  each  case  as  is  shown  in  sketch  B: 


Sketch  B 


This  fitting  although  made  visually  should  be  made  carefully.  If  the 
ordinate  of  pulse  1  is  designated  as  r^,  of  pulse  2  as  r^,  and  of 

pulse  n  as  rn,  then  the  time  response  to  a  unit  impulse  corresponding 
to  this  Re[H(io)j]  is  given  by 

h(t)  =  Yl  **  h  (t)  (6) 

j=l  J  J 


The  time  functions  h^t),  hg(t),  hn(t)  associated  with  each  of 

the  rectangular  pulses  shown  in  sketch  B  are  tabulated  in  reference  3. 
These  time  functions  are  the  inverse  Fourier  transforms  of  unit  rectan¬ 
gular  pulses  of  the  Re[H(iaO[]  . 

The  method  was  applied  to  the  frequency -response  data  of  figure  1. 
The  fitting  of  the  staircase  function  to  the  Re [H ( icu)]  is  illustrated 
in  the  upper  part  of  figure  4  and  the  accuracy  of  the  method  is  demon¬ 
strated  in  the  lower  part  of  figure  4.  Again  the  error  in  the  computed 
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response  to  the  unit  impulse  at  t  =  0  is  due  to  cutting  off  the 
Re[H(io))J  at  u>  =  20  radians/sec  and  not  integrating  to  °°  as  indi¬ 
cated  by  equation  (3). 


Fourier  Method 


Fourier  response  to  unit  impulse.-  In  reference  b,  the  response  to 
a  unit  step  input  is  derived  by  the  Fourier  method.  Differentiating 
this  response  gives  an  expression  by  which  the  response  to  a  unit  impulse 
may  be  approximated.  For  a  linear  system  characterized  by  its  frequency 


response  with  amplitude  ratio 


& 

8 


and  phase  angle  <f)^,  this  expression 


is 


h(t>=^r  If 

*  n=l  ° 


(2n-l)a>£ 


cos 


(2n  -  l)a*t  +  (Al^.D  J  (7) 


In  using  equation  (7)  the  choice  of  the  fundamental  frequency 


determines  the  accuracy  and  length  of  the  computations.  Since  the 
accuracy  is  affected  by  oip  it  has  been  found  by  experience  that  a 

—  — can  be  used  as  a  good  first  estimate,  where  un  i®  the 

5 


value  cup 


undamped  natural  frequency  of  the  system.  Instead  of  <%,  the  value  of 
frequency  at  which  the  amplitude  ratio  peaks  for  a  lightly  damped  system 
may  be  used.  This  value  may  be  determined  from  the  frequency  response 
of  the  system. 


In  table  I  the  numerical  computation  of  h(t)  for  the  system  defined 
by  the  frequency  response  of  figure  1  is  shown.  In  this  computation, 
l4  terms  were  carried  in  the  expression  for  h(t) ,  and  more  accuracy  may 
be  obtained,  of  course,  by  carrying  more  terms.  The  accuracy  of  the 
method  is  shown  in  figure  5  by  the  comparison  between  the  computed  and 
exact  response  to  the  unit  impulse. 


Other  Methods 

Schumacher's  method.-  The  method  of  reference  5  permits  the  computa¬ 
tion  of  transfer-function  coefficients  by  assuming  the  shape  of  the 
transfer  function  relating  the  input  and  output  and  then  curve  fitting 
this  relation  to  the  aircraft  data  in  frequency -response  form.  Once  the 
transfer-function  coefficients  are  known,  the  system  is  completely 
specified  since  the  response  to  a  unit  impulse  input  or  any  other  arbi¬ 
trary  input  may  then  be  computed  by  the  normal  methods  available  for 
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solving  differential  equations.  To  illustrate  the  method  a  sample 
computation  is  performed  by  using  the  system  defined  in  figure  1.  The 
differential  equation  relating  input  8  and  output  x  for  the  fre¬ 
quency  response  shown  in  figure  1  is 

D2x  +  KqDx  +  Kgx  =  K58  (8) 

By  applying  the  vector  least-squares  method  of  reference  5,  the  following 
set  of  simultaneous  equations  is  obtained: 


and 


B*  =  |!|sin  0x6 

These  equations  are  then  solved  simultaneously  for  the  transfer  coeffi¬ 
cients  Kj,  K-|_,  and  K2.  From  the  form  of  equation  (8),  it  is  seen  that 

the  response  to  a  unit  impulse  input  is  given  by  the  equation 


h(t) 


K3  -at  .  . 

=  —  e  sin  oat 

%  ^ 


(10) 
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where 


<m  = 


and 


a 


% 

2 


A  sample  computation  is  shown  in  table  II.  Table  II  demonstrates 
the  computational  steps  involved  in  solving  for  the  transfer  coeffi¬ 
cients  K2,  and  for  the  system  defined  in  figure  1.  These 

transfer  coefficients  are  then  substituted  into  equation  (10),  and  a 
plot  of  this  function  is  compared  with  the  exact  response  of  the  system 
to  a  unit  impulse  input  in  figure  6. 

This  method  gave  additional  information  when  applied  to  the  second- 
order  system  since  the  transfer -function  coefficients  were  determined  in 
the  process  of  the  computation.  A  reasonable  amount  of  computing  time 
was'  required,  and  good  accuracy  was  obtained.  The  method,  however, 
requires  previous  knowledge  of  the  form  of  the  transform  function 
relating  the  input  and  output.  Use  of  the  curve-fitting  method  on  the 
higher  order  system  did  not  prove  too  effective  in  determining  all  the 
parameters  of  the  flexible  system.  If  the  form  of  the  transfer  function 
of  the  short  period  is  assumed  and  fitted  to  that  portion  of  the  aircraft 
frequency  response  associated  with  the  short  period  then  the  method  is 
very  effective  in  determining  the  short-period  transfer-function 
coefficients. 

P-transform  method.-  The  P-transform  method,  as  described  in 
appendix  A,  differs  from  the  other  methods  presented  in  this  paper  in 
that  the  time  response  to  a  unit  impulse  can  be  determined  directly 
from  a  known  response  to  a  known  arbitrary  input.  The  computation  is 
carried  out  entirely  in  the  time  domain,  bypassing  the  frequency  plane 
entirely,  and  does  not  require  knowledge  of  the  transfer  function  relating 
the  input  and  output.  The  method  also  represents  a  simple  procedure  for 
using  the  response  to  a  unit  impulse  and  determining  the  response  to  a 
given  arbitrary  input. 

An  example  of  the  method  is  shown  in  figure  7  in  which  the  assumed 
input  S(t )  and  output  x(t)  are  given  from  which  the  response  to  the 
unit  impulse  h(t)  is  computed.  A  comparison  of  the  computed  response 
to  a  unit  impulse  and  the  exact  response  to  a  unit  impulse  is  shown  in 
the  lower  part  of  figure  7 • 
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Response  to  Arbitrary  Inputs 

The  procedure  for  obtaining  the  response  to  arbitrary  inputs  when 
the  response  to  the  unit  impulse  is  known  is  simply  a  matter  of  applying 
Duhamel's  integral  to  the  response  to  the  unit  impulse  and  time  history 
of  the  arbitrary  input.  .A  numerical  method  of  performing  this  operation 
is  shown  in  appendix  A  of  reference  3* 

The  response  to  a  unit  step  input  may  be  obtained  by  numerically 
integrating  the  response  to  a  unit  impulse  by  using  the  integrating 
matrix  given  in  reference  3.  In  like  manner,  the  response  to  a  ramp 
input  may  be  obtained  by  numerically  integrating  the  response  to  a  step 
function.  The'  response  to  a  triangular  input  is  obtained  by  superposi¬ 
tion  of  the  responses  to  various  ramp  inputs  translated  along  the  time 
scale. 

The  response  to  an  arbitrary  input  may  be  obtained  from  the  response 
to  a  unit  impulse  by  the  P-transform  method  by  use  of  equation  (Ak)  of 
the  appendix. 

For  directly  determining  the  time  response  to  an  arbitrary  input 
from  f r e quency-r e spons e  data,  the  Fourier  method  is  perhaps  the  best 
known.  The  general  method  is  indicated  in  references  6  and  7*  l11  general 

an  input  which  can  be  expressed  as  a  Fourier  series  can  be  represented  as 


00 

S(t)  =  A q  +  {si 

n=l  N 


sin  no^t 


(11) 


The  response  to  the  input  given  by  equation  (ll)  of  a  linear  system 
whose  amplitude  ratio  is  |~|  and  phase  angle  is  can  be  written  as 

00 

x  (t )  =  Bo  +  YZ  Cnmf 
n=l  x 

For  a  few  specific  shapes  of  inputs  the  values  of  B0,  Cna^,,  and 
of  equation  (12)  are  known.  Two  wave  forms  frequently  used  are  the 
square  wave  and  triangular  wave  and  examples  of  these  are  given. 


sin 


naif 


na¥t  +  ^  +  (0xs)n 


(12) 
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Response  to  a  square  wave.-  As  indicated  in  reference  4  the  Fourier 
series  for  a  unit  square-wave  input  5 (t)  may  be  written  as 


sin(2n  -  l)a3ft 


(13) 


and  the  response  to  this  unit  step  input  of  a  linear  system  of  which  the 
amplitude  ratio  is  |g-|  and  the  phase  angle  is  0xg  is  given  by 


sin 


(2n-l)a>ft  + 


(2n-l)o>f 


(1*0 


The  accuracy  of  the  computation  depends  on  the  choice  of  the  optimum 
fundamental  square-wave  frequency  and,  as  before,  a  suitable  value  is 

usually  odp  =  ^  where  o)q  is  the  lowest  undamped  natural  frequency  of 

5 

the  system. 

Response  to  a  triangular -wave  input . -  The  specific  form  of  the  input 
used  in  determining  the  time  response  to  a  triangular -wave  input  is  shown 
in  sketch  C: 


6(t) 


Sketch  C 

where  Tf  is  the  total  period  of  the  input,  Tj_  is  the  length  of  the 
base  in  seconds,  and  8  is  the  maximum  value  of  the  input. 

Without  giving  the  analytical  development,  representation  of  the 
Fourier  series  for  this  triangular -wave  input  is  given  by 
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00  [-  — 

S(t)  =2Z  C2n_1  sin  |(2n  -  l)o>ft  +  02n_1 


where 


^  =  radians/sec 
1 


(15) 


(16) 


c 


2n-l 


1  -  cos(2n  -  l)it 


(17) 


and 


^2n-l 


(2n  -  1) 


(18) 


* 


The  response  to  this  triangular -wave  input  of  a  linear  system  of  which 
the  amplitude  is  and  the  phase  angle  is  0  is  then 


x(t) 


uu 

-L 


n=l 


2n-l 


x 

6 


(2n-l)cof 


sin 


(2n  -  l)coft  + 


^(Sn-Dojf  +  (4S)  (2n.1)oi 


(19) 


» 


Here  again  a  suitable  first  approximation  is 


For  a  higher  order  system  such  as  a  linear  system  with  several 
structural  modes  (as  indicated  by  peaks  in  the  amplitude -ratio  curve), 
the  choice  of  a  fundamental  frequency  becomes  difficult  inasmuch  as 

odd  multiples  of  adf  must  give  the  natural  frequency  and  the 

frequency  of  each  of  the  higher  structural  modes. 


-ill 
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Jf  an  input  is  nonperiodic  and  cannot  be  expanded  into  a  Fourier 
series  directly,  the  Fourier  transform  of  the  input  is  used.  The  tables 
of  reference  3  may  be  used  to  determine  numerically  the  Fourier  transform 
of  an  arbitrary  input.  The  procedure  for  computing  the  response  of  a 
linear  system  to  such  an  arbitrary  input  is  briefly  described  as  follows. 
The  Fourier  transform  of  the  output  is  formed  by  multiplying  the  Fourier 
transform  of  the  arbitrary  input  by  the  frequency  response  of  the  system, 
which  is  the  ratio  of  the  Fourier  transforms  of  the  output  and  input. 

The  inverse  Fourier  transform  of  this  product  is  then  taken  which  is  the 
time  response  of  the  system  to  the  arbitrary  input.  The  tables  of  refer¬ 
ence  3  may  also  be  used  to  perform  the  operation  of  the  inverse  Fourier 
transform. 


COMPARISON  AND  USE  OF  METHODS 
Comparison  of  Methods 


A  comparison  of  the  methods  as  applied  to  the  second-order  system, 
defined  by  the  frequency  response  shown  in  figure  1,  may  be  made  by 
noting  the  differences  between  the  exact  response  to  a  unit  impulse  and 
the  response  computed  by  each  of  the  methods  shown  in  figures  2  to  7 • 
From  these  figures  there  appears  to  be  little  difference  between  the 
accuracy  of  any  one  method  over  the  others ,  and  in  each  case  greater 
accuracy  may  be  obtained  at  the  cost  of  more  computation.  For  the 
comparisons  shown,  however,  the  computing  time  required  by  the 
rectangular -pulse  method  was  significantly  less  than  the  time  required 
by  the  other  methods.  The  Fourier  response  to  a  unit  impulse  required 
the  most  time.  A  listing  of  the  methods  in  the  order  of  the  computing 
time  required  for  the  comparisons  shown  in  figures  2  to  7  is  given  as 
follows : 

Rectangular  pulse  (ref.  3) 

Schumacher ' s  (ref .  5 ) 

Floyd's  (ref.  2) 

Numerical  integration 

Fourier  response  to  a  unit  impulse 

(All  computations  were  performed  on  a  desk- type  computer.) 

A  more  severe  test  of  the  methods  occurs  when  they  are  applied  to 
linear  systems  of  higher  order.  In  order  to  demonstrate  an  application 
of  this  type,  a  further  comparison  of  the  methods  was  made  by  applying 
them  to  the  system  defined  by  the  transfer  function 


100 


225 


2  +  6s  +  10  s2  +  0.4s  +  100  s2  +  0.2s  +  225 


(20) 
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The  frequency  response  of  this  system  is  shown  in  figure  8.  The  real 
part  of  the  frequency  response  shown  in  figure  9  was  computed  from 
figure  8  by  using  the  relation 


x 

8 


cos  0x6 


(21) 


In  figure  10  a  comparison  of  the  response  to  a  unit  impulse  computed  by 
four  of  the  methods  is  shown.  For  Floyd's  method  a  27-line  fit  to  the 
Re  [~H(i  o))J  shown  in  figure  9  was  used.  For  the  rectangular -pulse  method, 

the  staircase  fit  to  the  Re[H(io))J  was  made  with  the  interval 
Aco  =  1  radian/sec,  and  for  the  numerical-integration  method,  an  interval 
Ao)  =  0.5  radian/sec  was  used.  The  response  of  the  linear  system  defined 
by  equation  (20)  to  a  triangular  input  shown  in  figure  7  of  reference  3  was 
used  to  compute  the  points  shown  in  figure  10  for  the  P-transform  method. 

A  time  interval  of  At  =  0.1  second  was  used.  The  Fourier  response  to  a 
unit-impulse  method  was  found  to  be  impractical  in  this  case  because  of  the 
large  amount  of  computation  required  to  obtain  any  accuracy.  The  accuracies 
of  each  of  the  four  methods  illustrated  appear  to  be  equivalent. 


Use  of  the  Response  to  a  Unit  Impulse  as  a 
Flight-Data-Analysis  Technique 


As  a  further  comparison  of  the  methods  applied  to  a  higher  order 
system,  a  typical  longitudinal  maneuver  (M  =  0.82)  for  a  flexible  swept- 
wing  airplane  has  been  analyzed.  In  this  case  the  output  was  the 
pitching -velocity  response  S  at  the  center  of  gravity  of  the  airplane 
and  the  input  was  the  elevator  angle  Se.  The  analysis  of  this  maneuver 
serves  to  demonstrate  the  use  of  the  method  of  obtaining  the  time  response 
of  linear  systems  to  a  unit  impulse  from  frequency -response  data  as  a 
flight -data -analysis  technique.  The  Fourier  integrals  of  the  output  0 
and  input  5e  were  evaluated  by  using  automatic  computing  equipment  and 
the  methods  of  integration  of  product  curves  were  used.  The  frequency 
response  obtained  by  dividing  the  Fourier  integral  of  the  output  by  the 

?0 

Fourier  integral  of  the  input  is  shown  in  figure  11.  - 


The  Re 


(ico) 


computed  from  the  frequency  response  of  figure  11  is  shown  in  figure  12 
In  figure  13  a  comparison  of  the  response  to  a  unit  impulse  computed  by 
three  of  the  methods  is  shown.  For  Floyd's  method  a  37-line  fit  to  the 


Re 


shown  in  figure  12  was  used  and 
values  of  time  to  define  adequately  its  shape. 


was  computed  at  enough 
For  the  rectangular- 


pulse  method,  the  staircase  fit  to  the  Re 


was  made  with  the 
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interval  As  =  1  radian/sec  and,  for  the  numerical-integration  method, 
an  interval  of  Ad  =0.5  radian/sec  was  used.  In  the  case  of  the 
numerical  method,  points  are  shown  in  figure  13  only  at  time  intervals 
of  0.1  second  in  order  to  compare  accuracies.  Several  attempts  were 
made  to  compute  the  response  to  the  unit  impulse  by  the  P-transform 
method  using  the  original  elevator  input  6e(t)  and  output  § (t )  shown 
in  figure  14.  The  computations  were  made  at  time  intervals  ranging  from 
0.02  second  ^  At  ^  0.1  second  but  did  not  yield  satisfactory  accuracy. 

It  appears  that  this  inaccuracy  is  primarily  due  to  the  sensitivity  of 
the  method  to  small  errors  in  the  first  few  terms  of  the  response  or 
input. 

Time  histories  of  the  original  elevator  input  6e  and  pitching- 
velocity  response  0  are  shown  in  figure  l4 .  The  elevator  input  5e 
has  also  been  used  as  a  forcing  function  with  the  response  to  a  unit 
impulse  shown  in  figure  13  to  compute  a  time  history  of  0.  The  Duhamel 
method  outlined  in  appendix  A  of  reference  3  with  a  time  interval  of 
0.05  second  was  used  for  these  calculations.  A  comparison  of  this 
computed  0  response  with  the  original  0  response  from  a  flight 
test  is  shown  in  the  lower  part  of  figure  14.  This  comparison  shows  the 
amount  of  error  involved  in  the  total  computation  procedure  (transferring 
the  data  in  time-history  form  to  frequency -response  form,  then  to  the 
response  to  a  unit  impulse,  and  finally  to  the  response  to  an  arbitrary 
input  by  means  of  Duhamel' s  integral). 

Another  flight-test  maneuver  at  M  =  0.80  was  selected  with  the 
other  conditions  approximately  the  same  as  the  previous  data  in  order 
to  see  how  well  the  response  to  a  unit  impulse  computed  from  one  maneuver 
could  be  used  to  predict  the  time  response  to  an  arbitrary  input  from  a 
different  maneuver.  The  time  histories  of  this  maneuver  are  shown  in  fig¬ 
ure  15.  Also  shown  in  the  lower  part  of  figure  15  is  the  0 (t )  response 
calculated  by  application  of  Duhamel 's  integral  to  the  response  to  a  unit 
impulse  given  in  figure  13,  and  the  elevator  motion  for  the  maneuver  is 
shown  in  the  upper  part  of  figure  15.  A  comparison  of  the  computed  and 
measured  0  response  in  the  lower  part  of  figure  15  gives  some  indication 
of  how  well  the  motions  of  an  aircraft  can  be  predicted  by  a  detailed 
analysis  of  a  single  maneuver. 


CONCLUDING  REMARKS 


Several  methods  were  compared  for  obtaining  the  time  response  of 
linear  systems  to  either  a  unit  impulse  or  arbitrary  input  from  frequency- 
response  data.  The  methods  were  compared  on  the  basis  of  accuracy, 
computing  time  required,  and  applicability  to  higher  order  linear  systems. 
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The  application  of  each  of  the  methods  to  a  simple  second-order 
system  indicated  little  difference  between  the  accuracy  of  one  method 
over  the  others,  and,  in  general,  it  would  be  expected  that  greater 
accuracy  could  be  obtained  for  each  of  the  methods  at  the  cost  of  more 
computing  time. 

For  higher  order  systems  the  three  methods  generally  classified  as 
inverse  Laplace  transform  methods  were  most  effective.  They  gave  good 
results  for  a  moderate  amount  of  computation.  All  of  these  methods  are 
based  on  the  evaluation  of  the  form  of  the  inverse  Laplace  transform 
equation  for  positive  values  of  time. 

The  methods  generally  classified  as  Fourier  methods  gave  good 
accuracy  when  applied  to  the  second-order  system.  The  accuracy  of  these 
methods  was  found  to  depend  on  the  choice  of  a  fundamental  frequency  a 

For  simple  systems  a  suitable  value  was  found  to  be  ouf  =  ^  where 

is  the  lowest  natural  frequency.  For  higher  order  systems  the  choice 
of  (Of  becomes  more  difficult.  The  Fourier  response  to  a  unit-impulse 
method  was  found  to  be  impractical  in  the  case  of  higher  order  systems 
because  of  the  large  amount  of  computation  required  to  obtain  accuracy. 

Schumacher’s  method  gave  additional  information  when  applied  to  the 
second-order  system  since  the  transfer -function  coefficients  were  deter¬ 
mined  in  the  process  of  the  computation.  This  method  required  a  reason¬ 
able  amount  of  computing  time  and  gave  good  accuracy.  The  method, 
however,  requires  previous  knowledge  of  the  form  of  the  transfer  func¬ 
tion  relating  the  input  and  output.  Use  of  the  method  on  higher  order 
systems  did  not  prove  too  effective  in  determing  all  the  parameters  of 
the  flexible  system.  A  valuable  use  of  the  method  was  found  in  fitting 
the  known  short -period  transfer  function  to  only  the  short-period  portion 
of  the  flexible -system  frequency  response  to  determine  the  short-period 
transfer  coefficients. 

The  P-transform  method  is  different  from  the  other  methods  presented 
in  that  the  time  response  to  a  unit  impulse  can  be  determined  directly 
from  a  known  response  to  a  known  arbitrary  input.  The  computation  is 
carried  out  entirely  in  the  time  domain  and  bypasses  the  frequency  plane 
entirely.  The  method  also  represents  a  simple  procedure  for  using  the 
response  to  a  unit  impulse  and  determining  the  response  to  a  given 
arbitrary  input.  The  method,  however,  when  applied  to  the  flight  data 
of  this  paper  did  not  yield  satisfactory  accuracy.  It  appears  this  is 
primarily  due  to  the  sensitivity  of  the  method  to  errors  in  the  first 
few  terms  of  the  response  or  input. 

For  all  the  examples  computed  in  this  paper  the  rectangular-pulse 
method  of  MCA  Technical  Note  3598  required  less  computing  time,  for  the 
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same  degree  of  accuracy,  than  any  of  the  other  methods  used  to  obtain 
the  time  response  to  a  unit  impulse  from  frequency-response  data. 

When  an  aircraft  frequency  response  is  known,  it  appears  that  these 
methods  may  be  used  as  a  flight-data-analysis  technique  which  permits 
prediction  of  aircraft  motions  and  loads  without  knowledge  of  the  equa¬ 
tions  of  motion  relating  the  inputs  and  outputs  for  a  flexible  aircraft. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  March  l6,  1958. 
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APPENDIX  A 

P -TRANSFORM  METHOD 


In  reference  8  a  linear  operational  calculus  is  introduced  which 
appears  to  be  well  adapted  to  the  numerical  analysis  and  synthesis  of 
linear  systems .  In  this  calculus  a  polynominal  transform  or  P-transf orm 
of  a  f inaction  F(t)  is  defined  by  the  equation 


P|j(t  jj  =  F(x) 

=  YZ  F(mAt)xm  (Al) 

m=-K 


where  F(mt)  are  the  ordinates  of  F(t)  at  integral  multiples  of  a 
time  interval  At.  The  inverse  P-transf  orm  is  given  by 


F(t)  =  P-1_F(x)] 


(A2) 


The  superposition  (convolution  or  Faltung)  integral,  which  is  also  known 
as  Duhamel's  integral,  relates  the  input  b(t),  output  x(t),  and 
response  to  a  unit  impulse  h(t)  (also  called  a  memory  function)  of  a 
linear  system  and  is  given  by 


x(t)  =  f  6(t)  h(t-r)dT  (A3) 


Bubb  (ref.  8)  shows  that  the  P-transf  orm  of  equation  (A3)  is 

x(x)  =  At  b(x)  h(x)  (A4) 


In  the  synthesis  problem  the  input  5(t)  and  the  output  x(t)  are 
given  and  the  response  to  the  unit  impulse  h(t)  is  to  be  calculated. 
This  calculation  is  performed  by  forming  the  P-transf orms  of  x(x)  and 
b(x)  and  dividing  by  ordinary  polynominal  division,  x(x)  by  S(x),  to 
get  the  P-transform  of  the  response  to  the  unit  impulse 
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At  5(x) 

L 


(A5) 


For  practical  cases,  however,  the  summation  ^is  made  over  a  finite  range 
of  t  values.  The  inverse  P-transform  to  h(x)  then  gives  h(t),  the 
time  response  to  the  unit  Impulse.  An  example  of  this  operation  is  shown 
in  table  III;  the  known  input  6(t)  and  output  x(t)  are  shown  in  the 
upper  part  of  figure  7,  and  the  computed  response  to  the  unit  impulse  is 
compared  with  the  analytic  solution  in  the  lower  part  of  figure  7. 

This  operational  calculus  is  also  well  adapted  to  the  solution  of 
the  analysis  problem  in  which  the  system  response  to  a  unit  impulse  h(t) 
and  the  input  function  &(t)  are  known  and  the  calculation  of  the  system 
output  x(t)  is  desired.  The  P-transforms  of  h(t)  and  6(t)  are 
formed  and  multiplied  together  by  ordinary  polynominal  multiplication  as 
indicated  in  equation  (A4).  Since  this  operation  is  just  the  inverse  of 
the  operation  shown  in  table  III,  an  illustrative  example  is  not  shown. 
This  method  has  been  found  to  be  a  simple  and  rapid  means  of  applying 
the  Duhamel  process . 

The  value  of  this  operational  calculus  lies  in  the  fact  that  all 
computations  remain  in  the  time  domain  and  no  translation  to  the  frequency 
plane  is  required.  Also,  only  simple  direct  arithmetical  procedures  are 
required  for  solving  practical  problems. 
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TABLE  I 

NUMERICAL  CALCULATION  ILLUSTRATING  COMPUTATION  OF  THE 
FUNCTION  FOR  THE  RESPONSE  TO  A  UNIT  IMPULSE 
BY  TEE  FOURIER  METHOD 


Oip  =  ^  =  0.64 


n 

2n  -  1 

(2n  -  l)cDf 

X 

8 

0x5  > 

deg 

1 

1 

0.64 

0.0970 

-22 

2 

3 

1.92 

.O76O 

-61 

3 

5 

3.20 

.0520 

-91 

4 

7 

4.48 

.0350 

-111 

5 

9 

5.76 

.0240 

-124 

6 

11 

7.04 

.0170 

-133 

7 

13 

8.32 

.0130 

-140 

8 

15 

9.60 

.0100 

-145 

9 

17 

10.89 

.0079 

-149 

10 

19 

12.16 

.0064 

-152 

11 

21 

13.45 

.0053 

-155 

12 

23 

14.72 

.0044 

-157 

13 

25 

16.00 

.0038 

-159 

14 

27 

17.50 

.0033 

-160 

h(t)  =  ——[0.O97  cos(0.64t  -  22)  +  O.O76  cos(l.92t  -  6l)  + 

jt  L 

0.052  cos(3.2t  -  91)  +  0.035  cos(4.48t  -  111)  + 

0.024  cos(5.76t  -  124)  +  0.017  cos(7.04t  -  133)  + 
0.013  cos(8.32t  -  140)  +  0.010  cos(9.6t  -  145)  + 
O.OO79  cos(l0.89t  -  149)  +  0.0064  cos(l2.l6t  -  152)  + 
0.0053  cos(l3.45t  -  155)  +  0.0044  cos(l4.72t  -  157)  + 
0.0038  cos(l6t  -  159)  +  0.0033  cos (17 .3t  -160)] 
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TABLE  II 

NUMERICAL  CALCULATIONS  ILLUSTRATING  SCHUMACHER'S  METHOD' 


radians/sec 

in 

0x5 } 
deg 

Ax  =  |||  cos  0x£. 

=  III  Sin  ^x5 

0 

0.1000 

0 

0.10000 

0 

1 

.0924 

-33.7 

.07687 

-.03127 

2 

.0745 

-63.4 

.03336 

-.06661 

3 

.0555 

-86.8 

.00310 

-.05541 

4 

.0404 

-104.0 

-.00977 

-.03920 

5 

.0298 

-116.6 

-.01334 

-.02664 

6 

.0225 

-125.8 

-.01316 

-.01825 

7 

.017^  . 

-132.9 

-.01184 

-.01274 

8 

.0139 

-138.4 

-.01039 

-.00923 

k3  =  1.0018 
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tICAL  CALCULATION  ILLUSTRATING  THE  P-TRANSFORM  METHOD 


2k 


tx  ?(D 


<D  O 
W  Q-i 
U  W 

II 


a3  ^ 
w  <D 

o  w 

0)  ft- 

M  g 

§5' 
ft  -P 

Cfl  -H 

&  § 


cd  * 

v 

O  tfi 

41 3 

its 

o 

ft  4? 

m 


<v  d 

«  3 


$■$ 

g* 


t*  'x 


— 

<3 * 


*§5 

a'tio 


.4-  CM  VO  ON  t—  VD  IfN  t—  VO  KN  . 
H  t~0  H  H  O  ONtr-VD  LTV. 
OOHHr-lHOOOO 


o  ir\  h 

fA  CVI  OJ 

o  o  o 


I 


OH  CM  KN  _rf  1A  VO  t-CO 
CM  KN  ,4  IT\  VO  t“  CO  ON  H  H  H  H  H  H  HHH 
XXXXXXXXXXXX  X  X  X  XXX 
wvoo\t^v£)inb“VOfn4-oinH  in  vd  vo  j  o 

~  ‘  CM  CM  H  H  O  O  KN 

o  o  o  o  o  o  o 


CM  VO  CTv  t —  vo 
H  t-O  H  H  O  ON  t-  VO  ITN 
OOHHHHOOOO 


O  H  CM  KN  ,4  irvvor-CO 
CM  NN  41  n  VO  tN-^OO  ONHH  H  h.h  h  hh  h 
XXXXXXXXXXXX  X  X  x 
OjfAHt-^vOVOnO  t~-  H  KN  H 


oovOfAfngVHvominNco  on 
HajvooinOVfAj-j'aj  On  vo  k\ 
O  O  O  H  H  H  CM  CM  CM  CM  H  r— !  H 


o 

o  _  _ 

o  o  o  o 


OOOOOOOO  o 


KN 

CO  CT\  KN 

rH 

o  o 


£vg  eg 

O  rl  lA 
ON  C-  £> 

8  0  o 
o  o 


S3 


KN  ITN  ^  KN 


OOOOOOOOO 


4  O 

H  O  ON  O  KN  KN 
OH  CM  VO  O  UN  ON  rA  -4  -4  CM  ON  VO 

8  000HHHCMCMCM0JH  H 

ooooooooooo  o 


KN  VO  CM 

o  h  in 

UN  C—  IfN 
O  O  O 
O  O  O 


H  CM  KN  4  iTn  4  KN  CM  H 


OOOOOOOOO 


H  OJ  lA  4  liAVO  t—  CO  ON  O  H  CM  KN  -4  ITN  VO  t—  CO 
O  HHH  H  HH  HHH 


NA.CA  TN  3701 


JX  fetO 


^X 

*3 


Response  to  a  unit  impulse,  h(t)  Re[H(icjJ| 


26 


MCA  TN  3701 


Figure  2.-  Time  response  to  a  unit  impulse  obtained  from  Floyd's  method 
compared  to  the  exact  values  for  the  second-order  system. 
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Figure  5.-  Time  response  to  a  unit  impulse  obtained  from  the  Fourier 
method  compared  to  the  exact  values  for  the  second-order  system. 


Figure  6.-  Time  response  to  a  unit  impulse  obtained  from  Schumacher* 
method  compared  to  the  exact  values  for  the  second-order  system. 
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Figure  9--  Real  part  of  the  frequency  response 

order  system. 
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Figure  11.-  A  typical  frequency  response  obtained  from  flight  data  of  a 

flexible  airplane  at  M  =  0.82. 
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for  the  typical  frequency  response  obtained  from  flight  data  of 
flexible  airplane  at  M  =  0.82. 
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Figure  14 . -  Comparison  with  the  actual  flight  record,  of  the  calculated 
pitching-velocity  response  using  the  same  input-output  data  as  was 
used  to  obtain  the  typical  frequency  response  of  the  flexible  air¬ 
plane  at  M  =  0.82. 
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predicting  loads  and  motions  for  a  flexible  aircraft  predicting  loads  and  motions  for  a  flexible  aircraft 

when  the  aircraft  frequency  response  is  known.  when  the  aircraft  frequency  response  is  known. 


